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Abstract 

This paper is concerned with decay estimate of solutions to the semilin- 
ear wave equation with strong damping in a bounded domain. Introducing 
an appropriate Lyapunov function, we prove that when the damping is 
linear, we can find initial data, for which the solution decays exponentially. 
This result improves an early one in [4]. 



1 Introduction 

In this paper we are concerned by the following problem 

Utt ~ Aw — ujAu t + (MUt = u\u\ p ~ 2 x € CI, t > 
u(x,t) =0, ie9Sl,i>0 (1.1) 

u(x, 0) = uq(x), ut(x,Q) = u\{x) x G 

in a bounded regular domain fl C K n . Here p > 2 and arc positive con- 
stants. Only one of this constant must be strictly positive (as pointed out 
by Gazzola and Squassina [4]). We will suppose that fj, > and u > (see 
Remark 2.4 for the case uj = 0). The present problem has been studied by Gaz- 
zola and Squassina [4]. In their work, the authors proved some results on the 
well-posedness and investigate the asymptotic behavior of solutions of problem 
(1.1). In particular, they showed the global existence and the polynomial decay 
property of solutions provided that the initial data are in the potential well, [4, 
Theorem 3.8]. The proof in [4] is based on a method used in [8] and [7]. In 
these works, the authors obtain the following differential inequality: 



d_ 



(l+t)E(t) <E(t) 
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where E is the energy of the solution. Unfortunately they obtain a decay rate 
which is not optimal. 

The nonlinear wave equations related to (1.1) has been investigated by many 
authors [1, 3, 6, 8, 14, 15, 20, 21]. 

In the absence of the nonlinear source term, it is well known that the presence 
of one damping term i.e. oj > or fi > ensures global existence and decay of 
solutions for arbitrary initial data (see [6, 11]). For u = /i = 0, the nonlinear 
term u \u\ p causes finite-time blow-up of solutions with negative initial energy 
(see [1, 10]). 

The interaction between the damping and the source terms was first con- 
sidered by Levine [12, 13]. He showed that solutions with negative initial en- 
ergy blows up in finite time. When uj — and the linear term u t is replaced 
by u t \u t \ m ~ 2 , Georgiev and Todorova [5] extended Levine's result to the case 
where m > 2. In their work, the authors introduced a method different from 
the one known as the concavity method. They determined suitable relations 
between m and p, for which there is global existence or alternatively finite time 
blow-up. Precisely, they showed that the solution continues to exist globally "in 
time" if m > p and blows up in finite time if p > m and the initial energy is 
sufficiently negative. Vitillaro [19] extended the results in [5] to situations where 
the damping is nonlinear and the solution has positive initial energy. Similar 
results have been also established by Todorova [17, 18], for different Cauchy 
problems. 

We recall here that the potential well method introduced by Payne and 
Sattinger [16] is also useful and widely used in the litterature to investigate 
the local existence, global existence and asymptotic behavior of the solutions to 
some problems related to problem (1.1) (see [18, 19, 3, 15, 7, 8, 20]). Introducing 
a strong damping term Au t makes the problem different from the one considered 
in [5]. For this reason less results are, at the present time, known for the wave 
equation with strong damping and many problems remain unsolved (see [4]). 

The purpose of this paper is to obtain a better decay estimate of solutions to 
the problem (1.1). More precisely we show that we can always find initial data 
in the stable set for which the solution of problem (1.1) decays exponentially. 
The key tool in the proof is an idea of Haraux and Zuazua [6] and Zuazua [21], 
which is based on the construction of a suitable Lyapunov function. This kind 
of Lyapunov function, which is a small perturbation of the energy, has been 
recently used by Benaissa and Messaoudi [2] to study the exponential decay if 
a weakly damped semilinear wave equations. 

2 Asymptotic stability 

In this section, we introduce and prove our main result. For this purpose let us 
introduce the definition of the solution of problem (1.1) given by Gazzola and 
Squassina in [4]. 
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Definition 2.1. For T > 0, we denote 

Y _ f ue c ([o,n-ffo 1 (fi))n c 1 ([o,T] ) L 3 (n))nc a ([o ) T],F- 1 (n)) 
T I « ( ei 2 ([0,T],L 2 (l])) 

Given uq G -//q(O) and ui G L 2 (il), a function u £ Yt is a local solution to 
(1.1), if u(0) = Uo, u*(0) = Ui and 

Uu4>dx + / Vu\7(pdx + uj / \7utV(j)dx + /i / u t (f)dx = I \u\ p ~ 2 ucj)dx, 
Jn Jn Jn Jq 

for any function (f) G Hq(Q) and a.e. t G [0,T] 

Let us first define the Sobolcv critical exponent p as: 

{2iV 
— — -, for u> > and N > 3 
2jV _ 2 and p = oo, if TV = 1,2 . 
, for uj = and N > 3 
iV-2 ' 

We first state a local existence theorem whose proof is written by Gazzola and 
Squassina, [4. Theorem 3.1]. 

Theorem 2.1. Assume 2 < p < p. Let uo G Hq(H.) and ui G L 2 (VL). Then 
there exist T > and a unique solution of (1.1) over [0,T] in the sense of 
definition 2.1. 

As in the work of Gazzola and Squassina, [4] , we define the global solutions 
and the blow up solutions. 

Definition 2.2. Let 2 < p < p , uq G Hq(H.) and u% G L 2 (Vl). We denote u the 
solution of (1.1). We define: 

Tmax = sup|r > , u = u(t) exists on [0, T] j 

Since the solution u eYt (the solution is "enough regular"), let us recall that if 
T m ax < oo, then 

^ 15m ||Vu|| 2 + \\uth = +°o • 

If T max < oo, we say that the solution of (1.1) blows up and that T max is the 
blow up time. 

If Tmax — oo, we say that the solution of (1.1) is global. 

In order to study the blow up phenomenon or the global existence of the 
solution of (1.1), we define the following functions: 

L(u(t)) = \\Vu(t)\\ 2 ~\\u(t)r p , (2.1) 

JWi)) = i||Vu(t)||5-i||ti(f)||* (2.2) 
2 P 
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and 

E(u(t)) = J(u(t)) + ±\\u t (t)\\2 (2-3) 
To have a lighter writing of /, J and E, wc will write : 

I{t) = I(u(t)) , J(t) = J(u(t)) and E{t) = E(u(t)) 

Let us remark that multiplying (1.1) by ut, integrating over f2 and using inte- 
gration by parts we obtain: 

^ = -u,\\Vu t \\l-fi\\u t \\l,Vt>0. (2.4) 

Thus the function E is decreasing along the trajectories. As in [16], the potential 
well depth is defined as: 

d = inf maxJ(Ait). (2.5) 
«ei?i(n)\{o} a>o 

We can now define the so called "Nehari manifold" as follows: 

Af={u e H>;(n)\{0}; l(t) = 0}. 
J\f separates the two unbounded sets: 

N + = {u e i?o(0); I(t) > 0} U {0} and M~ = {u G H^(Q,);I(t) < 0} . 
The stable set W and unstable set U are defined respectively as: 
W = {u e H%(Cl); J(t) < d} nAf+ and U = {u G H£(Q,); J(t) < d} n M~ . 
It is readily seen that the potential depth d is also characterized by 

d = min J (u) . 

As it was remarked by Gazzola and Squassina in [4] , this alternative character- 
ization of d shows that 



/? = dist(0,A0 = min||Vu||2 = 4/-^ >0 . (2.6) 

uGAf y p — 2 

In the lemma 2.1, we would like to prove the invariance of the set J\f + : if the 
initial data uo is in the set M + and if the initial energy E(0) is not large (we 
will precise exactly how large may be the initial energy), then u{t) stays in N + 
forever. 

For this purpose, as in [4, 19], we denote by C* the best constant in the 
Poincarc-Sobolev embedding Hq(VI) <-» L P (Q) defined by: 

C" 1 =inf{||Vu|| 2 :7ie J ff 1 (r!),||w|| p = l} . (2.7) 



Let us remark (as in [4, 19]) that if p < p the embedding is compact and the 
infimum in (2.7) (as well as in (2.5)) is attained. In such case (see, e.g. [16, 
Section 3]), any mountain pass solution of the stationary problem is a minimizcr 
for (2.7) and C* is related to its energy: 

d = V — C- 2p/{p - 2) . (2.8) 
2p 

Remark 2.1. It is well know from the potential well theory, [16, 9] , that 
for every solution of (1.1), given by Theorem 2.1, only one of the following 
assumption holds: 

i) if there exists some to > such that u(to) G W and E(to) < d, then 
Mt > t , u{t) G W and E(t) < d. 

ii) if there exists some t a > such that u(t a ) G hi and E(t ) < d, then Vi > 
t , u(t) G U and E(t) < d. 

in) Vi > 0, E{t) > d . 

We can now proceed in the global existence result investigation. For this 
sake, let us state two lemmas: these two results are stated in [4, Proof of 
Therorem 3.8] but are not detailed. For a better understanding of the results, 
we give a short proof of these two results. 

Lemma 2.1. Assume 2 < p < p. Let uq G J\f + and u\ G L 2 (f2). Moreover, 
assume that E(0) < d. Then for any < T < T max , u(t, .) G Af + for each 
t G [0,T). 

Remark 2.2. Let us remark, that if there exists to G [0,T) such that 

E(t Q ) < d 

the same result stays true. It is the reason why we choose to = 0. 

Moreover , one can easily see that, from (2.8), the condition E(0) < d is 
equivalent to the inequality: 

P-2 

cl {jh m ) <: (2 - 9) 

This last inequality will be used in the remaining proofs. 

Proof. Since I(uo) > 0, then by continuity, there exists T* < T such that 
I(u(t, .)) > 0, for all t G [0,T»). Since we have the relation: 

P-2 |ln , |2 , 1 



we easily obtain : 



2p p 



J(t)>P— 1| Vu||l, ViG[0,T» 
2p 
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Hence we have: 



p * 



From (2.2) and (2.3), we obvioulsy have Vi G [0,T*), J(t) < E(t). Thus we 
obtain: 



\Vu\\l < -^-E(t) 
p-2 



Since E is a decreasing function of t, we finally have: 

||v«||l<-%J?(o),vtG[o l r.). (2.10) 

By definition of C*, we have: 

||«||?<CJ||V«||S = CJ||Vu||r a ||V«||i 
Using the inequality (2.10), we deduce: 

P-2 

NIS < c? (^tf(o)J 2 ||v«||l . 

Now exploiting the inequality on the initial condition (2.9) we obtain: 

Nl?<liv«||l ■ (2.11) 

Hence || V«||i-||«||P > 0, Vi G [0, T»), this shows that u (t, .) GAT+, Vie [0,T*). 
By repeating this procedure, T* is extended to T . □ 

Lemma 2.2. Assume 2 < p < p. Let uq G W + and G L 2 (fl). Moreover, 
assume that E(0) < d. Then the solution of the problem (1.1) is global in time. 

Proof. Since the map t <— >■ E(t) is a decreasing function of time t, we have: 



E(0) > E{t) = \\\u t \\l + ?1||V U ||| + -I(t) , 
2 2p p 



which gives us: 



s(o)>i|K||5 + ^-^||v«||l 



2" 2p 

Thus, V< G [0, T) , || Vu||2 + ||ut||2 is uniformcly bounded by a constant depending 
only on E(0) and p. Then by definition 2.2, the solution is global, so T max = 
oo. □ 

We can now state the asymptotic behavior of the solution of (1.1). 

Theorem 2.2. Assume 2 < p < p. Let u G W+ and u\ G L 2 (fl). More- 
over, assume that E(0) < d. Then there exist two positive constants C and £ 
independent of t such that: 

< E(t) < Ce~ & , Vf > 0. 



G 



Remark 2.3. Let us remark that these inequalities imply that there exist positive 
constants K and £ independent of t such that: 

l|Vu(t)||! + |Mt)||2 <Ke<\ Vi>0. 

Thus, this result improves the decay rate of Gazzola and Squassina [4, Theorem 
3.8], in which the authors showed only the polynomial decay. Here we show that 
we can always find initial data satisfying uq G J\f + and U\ G L 2 (Vl) which verify 
the inequality (2.9), such that the solution can decay faster than l/t, in fact 
with an exponential rate. Also, the same situation happens in absence of strong 
damping (to — 0). 

Proof. Since we have proved that Vf > , u(t) G Af + , we already have: 

0<E(t) Vi>0. 

The proof of the other inequality relies on the construction of a Lyapunov func- 
tional by performing a suitable modification of the energy To this end, for 
e > 0, to be chosen later, we define 



L(t)=E(t) + e / u t udx+ — \\Vu\\ 



(2.12) 



It is straightforward to see that L(t) and E{t) are equivalent in the sense that 
there exist two positive constants (3\ and /3a > depending on e such that for 
t > 

P\E(t) < L(t) < p 2 E{t). (2.13) 

By taking the time derivative of the function L defined above in equation (2.12), 
using problem (1.1), and performing several integration by parts, we get: 



dL(t) 
dt 



Vt**||i — /*||t*t||l + — e|| V«||| 

-e||tt||^ — £|U / u t udx . 



(2.14) 



Now, wc estimate the last term in the right hand side of (2.14) as follows. 
By using Young's inequality, we obtain, for any 5 > 



1 

45 



Utudx < — \\ut\\% + 8\\u\\% 



(2.15) 



Consequently, inserting (2.15) into (2.14) and using inequality (2.11), we have: 

dL ^) / ,.„V7.. ,12 



dt 



< 



-w||Vut||£+(e 
( 

HC 2 J + CP 

V 



A5 



A*) \Wth 



\ 



2p 



(P-2) 



E{0) 



IV«||2 



(2.16) 



<o 
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From (2.9), we have 

p-2 

°* ((^2) S(0) ) 2 - 1 < ° • 
Now, let us choose 8 small enough such that: 

P-2 

From (2.16), we may find rj > 0, which depends only on 8, such that: 
dL(t) 



dt 



< -oj\\Vu t \\l + (e + l) - M ) Kill - eTj||Vtt||l 



Consequently, using the definition of the energy (2.3), for any positive constant 
M, we obtain: 

^ < -M £j B(i)+(s(^ + l + ^)- M )|H|l- W ||V^||l 

y-»?)l|Vu||2 • (2.17) 
Now, choosing M < 2r], and e small enough such that 

£ (t8 + 1 + y)-») <0 - 



inequality (2.17) becomes: 

dL(t) 
dt 



< -MeE{t), Vi > 0. 



On the other hand, by virtue of (2.13), setting £ = — Me/ 02, the last inequality 
becomes: 

^ <-££(*), Vt>0 . (2.18) 

Integrating the previous differential inequality (2.18) between and t gives the 
following estimate for the function L: 

L(t) < Ce'V , Vt > . 
Consequently, by using (2.13) once again, we conclude 

E{t) < Ce - ** , Vi > . 
This completes the proof. □ 



Remark 2.4. Note that we can obtain the same results as in Theorem 2.2 in 
the case w = 0, by taking the following Lyapunov function 



L(t) = E{t) +e u t udx. 
Ja 

Remark 2.5. It is clear that the following problem: 

Utt — div ( — — ] — uiAu t + fJ,Ut = u\u\ p ~ 2 i£fl, t > 

Vv 1 + |Vu| 2 / 

u(x,t) = 0, x G on, t > 

, u(x, 0) = uo(x), ut(x, 0) = ui(x) x G n 

could be treated with the same method and we obtained also an exponential decay 
of the solution if the initial condition is in the positive Nehari space and its 
energy is lower that the potential well depth. 
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